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We report two new first-order Raman modes in the spectra of few-layer M 0 S 2 at 286 cm -1 and 
471 cm' 1 for excitation energies above 2.4 eV. These modes appear only in few-layer M 0 S 2 ; therefore 
their absence provides an easy and accurate method to identify single-layer M 0 S 2 . We show that 
these modes are related to phonons that are not observed in the single layer due to their symmetry. 
Each of these phonons leads to several nearly degenerate phonons in few-layer samples. The nearly 
degenerate phonons in few-layer materials belong to two different symmetry representations, showing 
opposite behavior under inversion or horizontal reflection. As a result, Raman active phonons exist 
in few-layer materials that have nearly the same frequency as the symmetry forbidden phonon of 
the single layer. We provide here a general treatment of this effect for all few-layer two-dimensional 
crystal structures with an inversion center or a mirror plane parallel to the layers. We show that 
always nearly degenerate phonon modes of different symmetry must occur and, as a result, similar 
pseudo-activation effects can be excepted. 

PACS numbers: 02.20.-a, 78.30.-j, 78.67.-n 


I. INTRODUCTION 

Two-dimensional crystals have received a lot of atten¬ 
tion recently as they can have novel physical, chemi¬ 
cal and mechanical properties not found in their bulk 
counterparts. 1,2 Besides graphene, layered transition 
metal dichalcogenides and especially molybdenum disul¬ 
fide (M 0 S 2 ) show great potential for novel nanoelectronic 
and optoelectronic devices. 3- ' For characterizing layered 
materials, Raman and photoluminescence spectroscopy 
have been established as viable tools 8-16 . We report here 
two new resonant first-order Raman modes for few-layer 
MoS 2 . These modes appear in the Raman spectra for 
excitation energies above 2.4 eV. At 2.7 eV excitation 
energy, they are distinctly above the noise level; their 
intensity is comparable to second-order Raman peaks 
[see Fig.1(a)]. At UV excitation, their intensity becomes 
even larger than the intensity of the typically examined 
A! x /A\ g and E'/E g modes [see Fig.1(b)]. Therefore they 
can potentially be used as a quick and reliable way to dis¬ 
criminate between single- and few-layer M 0 S 2 . Recently, 
similar Raman modes, which are either Raman inactive 
or forbidden in backscattering geometry in the bulk mate¬ 
rial, were reported in few-layer WSe 2 , TaSe 2 , and MoTe 2 , 
however a systematic derivation of all vibrational modes 
in few-layer crystals, including their symmetry and dis¬ 
placement patterns is still lacking. 17-20 In this paper we 
present a detailed derivation of these modes using group 
theory. As the results follow from symmetry considera¬ 
tions alone, such Raman modes are not a unique feature 
of transition metal dichalcogenides, but are generally ex¬ 
pected in all IV-layer systems with an inversion symmetry 
center or a horizontal mirror plane. 

All of the new Raman modes have in common that 
they are related to single-layer vibrations with the same 
frequency which are Raman inactive or are not observ¬ 
able in the specific scattering geometry. Two effects are 
responsible for their appearance in the Raman spectra 
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FIG. 1: Raman spectra of single (1L), bi- (2L) and 
trilayer (3L) M 0 S 2 on Si/Si02 excited with (a) 2.7 eV 
and (b) 3.8 eV. The new modes are the E' and E g 
modes around 286 cm -1 . Furthermore, an additional 
A Xg /A x mode appears in few-layer M 0 S 2 471 cm -1 , 
which is related to the A% single-layer mode. 


of few layer systems: (i) The point group of the few- 
layer material changes compared to the single layer and 
depends on the stacking order and layer number. There¬ 
fore, also the representations of the vibration patterns 
in the few-layer system and their corresponding Raman 
tensors can change, (ii) Each normal mode of the single¬ 
layer system leads to N normal modes in the iV-layer 
system. These modes are nearly degenerate in frequency 
in case of weak interlayer coupling. For a system with 
an inversion center or a horizontal mirror plane paral- 
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lei to the layers, they belong to two different represen¬ 
tations of the few-layer symmetry group, showing op¬ 
posite character under inversion or horizontal reflection. 
For instance, a single-layer normal mode which is odd 
(■ ungerade ) under inversion or horizontal reflection leads 
to both, even (g evade) and odd normal modes in the 
few-layer system. The even modes can now be Raman 
active, thus we observe a pseudo-activation of the inac¬ 
tive single-layer mode. The newly observed modes are 
thus not Raman inactive modes but regular Raman ac¬ 
tive modes of the few-layer material with nearly the same 
frequency as those forbidden in the single layer. 

To derive the vibration patterns of the few-layer nor¬ 
mal modes, we discuss the applicability of a simple lin¬ 
ear chain model combined with symmetry considerations 
and compare it to DFT calculations. The splitting of 
the single-layer modes into two different symmetries in 
the few-layer system will affect, besides Raman scatter¬ 
ing, also all other physical processes where the symmetry 
of the normal modes is critical, such as electron-phonon 
scattering. Furthermore it is reasonable to assume to 
find a correspondingly splitting into nearly degenerate 
states of different symmetry also for the electronic wave 
functions. 


II. VIBRATIONAL PROPERTIES AND RAMAN 
SCATTERING SELECTION RULES OF 
MOLYBDENUM DISULFIDE 

We will first discuss the vibrational properties and Ra¬ 
man selection rules for single-, bilayer and bulk (2if) 
M 0 S 2 , as well as the Raman activation of the two new 
Raman modes for the bilayer, before we treat the case of 
a general IV-layer system. 

A. Bulk molybdenum disulfide 

Bulk M 0 S 2 is a layered crystal, it consists of covalently 
bound single layers of M 0 S 2 , each formed by a layer of 
molybdenum atoms covalently bound to two layers of sul¬ 
fur atoms. The individual M 0 S 2 layers in the bulk crystal 
are bound by van-der-Waals interaction. The structure 
of the most common polytype 2 H consist of two paral¬ 
lel layers rotated relatively to each other by 7r/6, such 
that the sulfur atoms of one layer are directly below the 
molybdenum atom of the other layer. 2 .H-M 0 S 2 belongs 
to the Dsh symmetry group, with the c-axis perpendic¬ 
ular to the layers as the main rotational axis. The six 
atoms in the unit cell lead to 18 normal vibrations, which 
decompose at the T point into the following irreducible 
representations: 21 ’ 22 

T 2H = A lg ®2A 2u ®2B 2g ®B lu ®E lg ®2E lu ®2E 2 g®E 2u 

(!) 

The in-plane vibrations have E symmetries and are 
twofold degenerate, while the out-of-plane vibrations 


have A and B symmetries. The E 2g , E\ g and A\ g 
symmetries correspond to Raman active modes. To 
describe the scattering geometries we use the notation 
(ej, e s ), where (e s ) is the polarization vector of the 
incoming (scattered) light. We define a cartesian co¬ 
ordinate system with the x and y vectors parallel to 
layer planes and the z vector parallel to the c-axis. The 
A\ g mode can be observed in (x,x), ( y,y ) and (z,z) 
scattering geometries, while the E 2g mode can be ob¬ 
served in ( x,x ), ( y,y ) and ( x,y ) scattering geometries. 
In contrast, the Ei g mode requires a scattering geom¬ 
etry involving a z-component, i.e. ( x,z ), (y,z) and 
(z,z). 23 All first-order Raman modes in bulk have been 
determined experimentally by Raman spectroscopy: 24,25 
oj(E* g ) = 32 cm' 1 , u{E\ g ) = 287 cm" 1 , u{E^ g ) = 
383 cm -1 and u>{A\ g ) = 409 cm -1 . Furthermore, the 
two IR-active optical phonons were measured by reflec¬ 
tion spectroscopy: 24 co(Ei u ) = 384 cm -1 and co(A 2u ) = 
470 cm” 1 . The E\ g and E\ u modes are nearly degenerate 
due to the low interaction between the layers. 21 

B. Single-layer molybdenum disulfide 

Single-layer molybdenum disulfide shows D^h symme¬ 
try. The unit cell consists of one metal and two sulfur 
atoms. The representation of the nine normal modes at 
the T point for the single layer can be decomposed into 
the following irreducible representations: 26 

r“ oS2 = A[ © E" © 2 A 2 © 2 E' (2) 

Figure 2 (a)-(f) shows the vibration patterns of the nor¬ 
mal modes. Due to the fact that the layers in the bulk 
are only weakly coupled, we find for all normal modes of 
the single layer two phonon modes of the bulk; Davydov 
pairs with only slightly modified frequencies and same 
displacement pattern (one with zero phase shift between 
the layers and one with a phase shift by 7r, see below). 15,16 
The A \, E' and E” symmetries correspond to Raman ac¬ 
tive modes. The scattering geometries of these modes are 
equal to the scattering geometries of the corresponding 
Raman active modes A\ g , E 2g and E\ g in the bulk. The 
E" mode requires therefore a scattering geometry with a 
z-component to be observable. However as single-layer 
M 0 S 2 is a two dimensional crystal, it is difficult to realize 
such a scattering geometry experimentally in backscat- 
tering. Due to their similarity, it has become customary 
to associate the the A\ and E' optical phonon modes of 
the single layer with the A\ g and E 2g modes of the bulk. 

C. Bilayer molybdenum disulfide 

The unit cell of bilayer molybdenum disulfide is the 
same as the one of 2 H stacked bulk MoS 2 - As a result, 
the vibration patterns of the bulk unit cell and the bilayer 
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FIG. 2: (a)-(f) Normal mode vibration patterns and symmetries of single-layer M 0 S 2 (upper panels). The 
frequencies denote the typical range found in few-layer M 0 S 2 . (g)-(l) The corresponding in-phase and anti-phase 
combination normal modes for bilayer (bulk) M 0 S 2 (lower panels). 



FIG. 3: Correlation diagram of the symmetries of the 
normal modes for bulk (2 H), single- and bilayer MoS 2 - 


are equal [see Fig. 2 (g)-(l)]. However, the point groups 
are different, as the main rotation axis in bilayer is not 
six-fold like in bulk, but three-fold. The six-fold rotation 
symmetry of bulk M 0 S 2 results from a screw axis trans¬ 
formation, including the translation along the c-axis. As 
this transformation is not possible for a finite number 
of layers, bilayer M 0 S 2 belongs to the D 3 d point group. 
Due to the reduced symmetry, there are only four differ¬ 
ent representations in the bilayer instead of eight for the 
bulk. Figure 3 shows the correlation of the symmetries 
of the normal modes between bulk, single- and bilayer 
M 0 S 2 . The 18 normal vibrations decompose at the T 
point into the following irreducible representations: 26 


r“ oS2 = 3A lff © 3 A 2u © 3 E g © 3 E u (3) 

Alternatively, the normal modes of the bilayer can be 
easily constructed from the single-layer normal modes: 
For each single-layer normal mode there are two normal 
modes of the bilayer obtained by combining the single¬ 
layer vibration in-phase and anti-phase (phase shift by 
7r), as shown in Fig. 2. The in-phase and anti-phase com¬ 
binations always show opposite transformation behavior 
under spatial inversion and thus belong to two differ¬ 
ent symmetry representations, i.e. gerade and ungerade. 
On the other hand, they have nearly the same phonon en¬ 
ergy in case of weak interlayer coupling. This leads to the 
generation of a Raman active mode in the bilayer from 
a previously Raman inactive single layer mode: From 


the single-layer A 2 mode an IR-active A 2u mode and a 
Raman active A\ g mode are constructed in the bilayer 
[Fig. 2 (1)]. This has been experimentally observed for 
WSe 2 , TaSe 2 , and MoTe 2 in Ref. 17-20. As the A 2u bulk 
phonon in M 0 S 2 has been measured by IR spectroscopy 
and calculations of the single-layer A'.j phonon energy cor¬ 
respond very well to one of the newly observed modes in 
the bilayer (Fig. 1), we assign this new mode at 471 cm -1 
to the bilayer A\ g vibration. In the bulk material the 
corresponding vibration pattern has B 2g symmetry and 
is therefore Raman inactive. The A 2u phonon of the bi¬ 
layer is in bulk still described by an A 2u representation 
and remains IR active. 

For the E" single-layer mode we observe a similar ef¬ 
fect regarding the Raman scattering geometry. The E" 
mode is Raman active in the single-layer, however it can 
not be observed in backscattering geometry. In bilayer 
with D 3 d symmetry, we find now an E g and an E u mode 
originating from the single-layer E" mode [Fig. 2 (i)[. 
The E g mode of the bilayer is - like the E" mode of the 
single layer - Raman active, but in contrast also observ¬ 
able in backscattering. Again, we find a good agreement 
between the newly observed mode at 286 cm -1 (Fig.l) 
and calculations, as well as with Raman measurements 
of the bulk E\ g phonon. Therefore, we assign this new 
mode to the bilayer E g symmetry. 


III. NORMAL VIBRATION MODES AND 
THEIR SYMMETRY REPRESENTATIONS IN 
AN A-LAYER SYSTEM 

The concept of combining the single-layer normal 
modes in-phase and anti-phase can be extended to a Al¬ 
layer system with an inversion center or a mirror plane 
parallel to the layers. For the sake of comprehensibility, 
we neglect in the following the degeneracy of the .E-type 
normal modes, as it will not affect the results. Gen¬ 
erally, the number of normal modes is proportional to 
the number N of layers, as the number of atoms in the 
unit cell of a few-layer material increases linearly with 
N. For normal modes all atoms move periodically with 
the same frequency. Thus the few-layer normal modes 
can be constructed only by a superposition of identical 
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non layer-order changing: E, C n ,a v , a d 
amplitude vector n unaltered 



(b) layer-order changing: /, S„,C' 2/ new amplitude vector n'= In 
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TABLE I: Multiplication table (direct product) of the 9 
point groups that contain only non layer-order changing 
operations with C, and C\h■ All 12 point groups that 
contain an inversion center or horizontal reflection plane 
symmetry can be written as such a product. 
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FIG. 4: Illustration of the effect of (a) the non 
layer-order changing operations and (b) the layer-order 
changing operations on a 4-layer structure. 


single-layer vibration patterns with possibly different am¬ 
plitudes in each layer (assuming mode mixing effects to 
be negligible, due to weak interlayer interaction). These 
amplitudes form a vector n € St N , with the amplitude 
of the single-layer normal mode in the jth-layer as the 
jth-component rij of the vector n. The normal modes 
of the few-layer system, originating from a given single¬ 
layer normal mode, are now described by a set of vectors 
n b) g j TV}, which form a basis for the $t N , 

as the normal modes are required to be linearly indepen¬ 
dent. To find the symmetries of the few-layer normal 
modes and the number of modes with a certain sym¬ 
metry, it is not necessary to determine the n <l ' > (which 
requires to find and solve the dynamical matrix); instead 
it is sufficient to consider the general behavior of the 
n/d under the symmetry transformation of the few-layer 
point group, as shown in the following. 

Layered materials can not belong to the polyhedreal 
point groups, i.e., T,Td,Th,0, and Oh in Schoenflies no¬ 
tation, as those groups contain rotations, whose rotation 
axes are not parallel or perpendicular to the layer planes 
and would therefore not transform them into themselves. 
All layered materials must therefore belong to one of the 
remaining 27 point groups. The symmetry transforma¬ 
tions found in these groups can be divided into two sets, 
depending on their effect on the layer order in a few- 
layer structure: one set of symmetry operations that do 
not change the order of the layers, i.e., {E,C n ,a v ,ad} 
in Schoenflies notation, and one set with the symmetry 
operations that do change the layer order {i, S n , C 2 , ah} 
(see. Fig. 4). Among the 27 point groups there are 9 
groups that contain only non-layer-order changing oper¬ 
ations, i.e., C\,C n and C nv with n = 2, 3,4,6. Further¬ 
more, there are 12 groups that have an inversion center or 
a mirror plane parallel to the layers, thus being the point 
groups relevant for the few-layer systems considered here. 
These 12 groups can always be written as a direct prod¬ 
uct of one of the 9 groups without layer-order changing 
operations with Ci or C\h (see Tab. I). As a result the 
few-layer point groups have always a uniquely defined 
subgroup that contains only non-layer-order changing op¬ 


erations. In the following we will denote this group as the 
non-layer-order changing subgroup. 

Obviously, under all symmetry operations that do not 
alter the order of the layers, the TV-layer normal modes 
will transform like the corresponding single-layer vibra¬ 
tion pattern. For the layer-order changing operations, we 
can focus on the i and ah operations, as all further layer- 
order-changing operations can be written as a product 
of an element of the non-layer-order-changing subgroup 
with i or ah ( e.g. C 2 = a v ah )• As a result, the behav¬ 
ior under inversion or horizontal reflection will determine 
the behavior of all other layer-order changing operations. 

In the following, we will focus on TV-layer systems with 
inversion symmetry. All results can be directly trans¬ 
ferred to systems with horizontal mirror symmetry sim¬ 
ply by interchanging i and ah- We define now an operator 
E that inverts the order of the components of the vector 
it is applied to: 

Eu = E (vi, ...,Vj-i,Vj) = (vj,Vj- 1 , ...,ul) (4) 

Inverting the few-layer vibration pattern leads effec¬ 
tively to an inversion of the order of the amplitudes - 
just like the effect of E. However the inversion can also 
change the amplitude signs depending on the definition 
of the basis displacement patterns that are used to con¬ 
struct the displacement patterns from a vector n^K This 
effect is illustrated for a bilayer system in Fig. 5. 

Nevertheless, for an even layer number, it is always 
possible to find a basis where applying E to a vector will 
result in a vector describing the inverted displacement 
pattern, e.g. like in Fig. 5(b), by choosing all pairs of 
basis vectors that are transformed onto each other un¬ 
der inversion accordingly. In this case the vectors n M, 
will be eigenvectors of E with eigenvalue k that can only 
become ±1. The sign of k decides now how the cor¬ 
responding few-layer vibration pattern transforms under 
inversion. Obviously, for even layer numbers, the few- 
layer normal mode will be even ( gerade ) for positive k 
and odd ( ungerade ) for negative. 

For odd layer numbers we have to consider that the 
central layer will be transformed onto itself and not onto 
another layer. As a result the few-layer vibration pat¬ 
terns can have only the same transformation behavior un¬ 
der inversion as the single layer vibration pattern, when 
the central layer has an displacement amplitude unequal 
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FIG. 5: Schematic illustration of the effect of the 
choice of the basis displacement patterns. Two different 
bases are defined by the displacement patterns B\ and 
B 2 for (a) and B[ and B' 2 for (b). In each basis the 
same two displacement patterns (Pi and P 2 ) are 
constructed, due to the different bases the 
corresponding vectors n and n' are unequal. Applying 
£ to the vectors n and n! results for (a) in displacement 
patterns that are reflected along a horizontal mirror 
plane, or are inverted at an inversion center for (b). 


zero. Furthermore we have to use — £ instead of £ to 
find the vector that describes the inverted displacement 
pattern when the corresponding single-layer vibration is 
odd. On the other hand, a few layer vibration, where the 
central layer shows non-zero displacement can be only 
described by an vector n 1 -' 1 ' 1 with k = +1, as the ampli¬ 
tude of the central layer will be left unaltered by £. All 
remaining basis vectors can be chosen such that their be¬ 
havior under inversion matches the behavior of the cen¬ 
tral layer. As a result using such a basis we find that for 
k = +1, the few-layer vibration will now behave exactly 
as the single-layer vibration under inversion, whereas it 
shows opposite behavior for k = — 1. 

We now turn to the question how many of the N nor¬ 
mal modes created by a given single-layer normal mode 
have k = +1 or k= —1. This question can be answered 
independently of the inter- and intralayer coupling, as 
the eigenvectors of £ divide ^R N into two orthogonal sub¬ 
spaces depending on the eigenvalue. 27 As the vectors nW 
form a basis for their respective subspaces and simulta¬ 
neously form a basis for the absolute number of the 
normal modes with the same eigenvalue can not change, 
regardless the amount of coupling. Therefore we always 
find (A"| (ceiling function of y) normal vibrations with 
k = +1 and [(jJ (floor function of with n = —1. 

The symmetry of the few-layer normal modes can now 
be easily found. As the few-layer normal mode trans¬ 
forms under the non-layer-order changing transforma¬ 
tions exactly as the single-layer normal mode used for 
the construction of the few-layer normal mode, we deter¬ 
mine the reduced representation of the single-layer nor¬ 
mal mode in the non-layer-order changing subgroup of 
the few-layer point group. Starting from this represen- 


Gi - Gui x C, 

^FL (gerade) 

D 

n FL (ungerade) 

FIG. 6: Correlation diagram of the normal modes in 
the single-layer point group ( Csl ), the few-layer point 
group ( Cfl ), and the non-layer-order changing 
subgroup ( C su b )• Several normal modes of the single 
layer can reduce to the same representation of C su b, 
indicated by the dotted lines. In contrast always two 
normal modes of the few-layer system, one gerade and 
one ungerade , reduce to the same representation of 
Cfl- 

tation, the few-layer point group, which is the product 
group of the the non-layer-order changing subgroup with 
Ci, contains two representations, one gerade and one 
ungerade. These are the possible symmetries of the few- 
layer normal modes, see Fig. 6. 


A. Application to few-layer molybdenum disulfide 

The point group of few-layer M 0 S 2 with 2iZ-like stack¬ 
ing depends on the layer number, see above and Ref. 
22. Single-layer M 0 S 2 has a horizontal mirror plane; in 
odd layer numbers this symmetry is preserved, as a re¬ 
sult they belong to the point group. Structures with 
an even layer number have an inversion center instead 
of the horizontal mirror plane; they belong to the D^d 
point group. Both groups have the same non-layer-order 
changing subgroup C% v (see Tab. I). 

At the T point of single-layer M 0 S 2 we find six different 
phonons which transform according to the A\ , A 2 , E' and 
E" representations of the D%h point group. The reduced 
representation of these phonons in C? yu are E for E' and 
E" and A\ for A\ and A". To determine the possible rep¬ 
resentations of the phonons in a few-layer structure with 
even layer number, we search the representations in D 3 d 
that reduce to E and A\ in C^ v . These are E g and E u for 
E , and A\ g and A 2u for A±. For odd layer numbers, using 
the same approach, we search the representations in D%h 
that reduce to E and A± in C^y. These are E' and E" 
for E, and A\ and A 2 for A-\. The results of the previous 
section, concerning the number of normal modes with 
k = +1 and k = —1, allow us now to determine number 
of the phonons with a specific symmetry in the few-layer 
crystals that are constructed from each single-layer nor¬ 
mal mode (see Table II). Due to interlayer coupling the 
degeneracy of the few-layer phonons will be lifted. How¬ 
ever, for the iV-layer normal modes originating from the 
optical single-layer phonons this effect is quite small. For 
the TV-layer normal modes originating from the acoustic 
phonons, on the other hand, the interlayer coupling leads 
for all but one modes to a non-zero phonon energy. Due 
to their vibration patterns, these optical modes are de- 
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TABLE II: Number of T-point phonon modes and their 
representation of the point group of single-, few-layer 
and bulk M 0 S 2 . ojq gives the typical range of the 
phonon frequencies. Bold letters denote the symmetries 
of Raman-active phonons which are allowed in 
backscattering geometry with the incoming light 
propagating perpendicular to the layer plane. Stars (*) 
denote the newly observed Raman modes. 
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E" 

+1 

N 

2 
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- 

-1 

N 

2 

• E u 
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N 

2 

• E u 
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+1 
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2 
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•Ai' 

B 2 g 


- 

-1 

N 

2 
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N—l 

2 
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noted as shear and layer breathing-like modes and have 
been extensively reported in the literature . 26,28-30 

From the E" and phonons of the single-layer the 
general analysis discussed above naturally predicts nor¬ 
mal modes with E g /E' symmetry at 286 cm -1 (see Fig. 
2 (i) and Fig. 7) and A\ g /A ( symmetry at 471 cm -1 
[see Fig. 2 (1)] in few-layer M 0 S 2 . These phonons are 
Raman active and explain the the newly observed Ra¬ 
man modes in Fig. 1. From a symmetry point-of-view, 
these modes always appear in few-layer structures with 
TV > 1. Many of the transition metal dichalcogenides 
have a structure and stacking order like 2H-MoS-2- While 
the phonon energies will certainly change, their symme¬ 
tries and therefore also their Raman selection rules fol¬ 
low the same structure as discussed above. Thus similar 
modes to the newly observed few-layer modes are ex¬ 
pected and have partially been observed in some isostruc- 
tural materials . 1 ' -20 


B. Application to AB-stacked graphene 

Monolayer graphene belongs to the D§h point group, 
and has therefore an inversion center and a horizontal 
reflection plane. In few-layer graphene (with AB stack¬ 
ing), the six-fold rotation axis is lost; instead it shows 
a three-fold main rotation axis. Also the inversion cen¬ 
ter and horizontal reflection plane are not always pre¬ 
served in the few-layers: Few-layer graphene with an 
even layer number has only an inversion center. With an 


TABLE III: Number of T-point phonon modes of 
single-, few-layer graphene and graphite and their 
representation of the point group. u >0 is the typical 
range of the phonon frequencies. Bold letters denote the 
symmetries of Raman-active phonons which are allowed 
in backscattering geometry with the incoming light 
propagating perpendicular to the layer plane. 
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odd number of layers it has only the horizontal reflection 
plane. Therefore, the point group is £ 3 ^ for even layer 
numbers and £ 3 /, for odd layer numbers with TV > 1, 
like in M 0 S 2. 31 The subgroup that does not contain the 
layer-order-changing operations is in both cases C% v . The 
monolayer has six T point phonons of A 2 U , £23, Eiu and 
£23 symmetry. The reduced representation of the A 2 U 
and £23 modes in C 3 V is A 1 ; E\ u and £23 reduce to 
£. This leads to the same possible symmetries for the 
few-layer graphene phonons as for few-layer M 0 S 2 with 
the same point group (see Tab. III). Again, the acoustic 
phonons of single-layer graphene lead to the optical shear 
and layer breathing like modes . 32-35 The E 2g phonon, 
which is responsible for the G mode of graphene, has Ra¬ 
man active counterparts of E g /E' symmetry in the few- 
layer structures. Notably, the £23 mode leads to normal 
modes of A\ g jA! x symmetry in few-layer graphene that 
may be observable in Raman spectroscopy from a sym¬ 
metry point of view. However, to our best knowledge, 
no experimental observation of these modes has been re¬ 
ported so far. 


C. Constructing the TV-layer normal modes 

We will now show how to use a simple linear-chain 
model to construct the displacement patterns of the few- 
layer vibrational modes. The application of this model 36 
is well established for the low frequency modes in few- 
layer graphene and other layered materials . 26,30,33 In this 
model each layer is represented by a classical harmonic 
oscillator which interacts only with it nearest neighbors 
directly. Such a system is similar to a resonator with open 
ends, whose solutions are standing waves . 37 For a finite 
number of oscillators, those standing waves form an enve- 
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FIG. 7: The normal modes of five-layer M 0 S 2 
originating from the optical E" phonon of the single 
layer. They are constructed using the eigenvectors cM 
of the solution of the linear chain model. The dotted 
(dashed) box shows the combination of single-layer 
vibrations leading to lower (higher) interlayer coupling 
energy. The vibration patterns calculated by DFT are 
the same despite small mode mixing effects with modes 
of the same symmetry. 


lope for the eigenvectors c-'d of the solutions - indexed by 
i - whose components Cj describe the amplitude of the 
j -th layer. Furthermore the finite number of oscillators 
limits the minimum wavelength of the standing waves; as 
a result there are N solutions for an TV-layer system. 

The normal modes originating from the single-layer 
acoustic modes can be easily constructed from the vec¬ 
tors c^), as the for the single-layer acoustic modes all 
atoms show the same displacement. Thus the product of 
the single-layer acoustic mode displacement vector with 
the vector component c ^ can be used directly as the 
displacement vector for all atoms in the jth-layer. 

For the few-layer normal modes generated from op¬ 
tical phonons, the individual atoms in each layer show 
different displacement vectors. Therefore it is necessary 
to consider which combination of the single layer vibra¬ 
tion pattern leads to a higher or lower interlayer coupling 
energy for neighboring layers. This can be easily de¬ 
rived from a bilayer system: Starting from the top layer 
(j = 1) we construct the vibration pattern of the few- 
layer system step by step, by applying the single-layer 
vibration pattern to the next layer such that it minimizes 
(maximizes) the interlayer coupling when the sign of the 
components c) of the neighboring layers are the same 
(opposite). In Figure 7 the construction of the normal 
modes originating from the E" (286 cm -1 ) single-layer 
phonon leading to the newly observed Raman modes is 
demonstrated for five-layer M 0 S 2 . The symmetries of the 
few-layer systems can be easily attributed to the two pos¬ 
sible representations by optical inspection, see above. As 
expected, we find two Raman active E' modes and three 
inactive E" modes. 

We have calculated the normal modes of the layered 
transition metal dichalcogenide WS 2 by density func¬ 
tional theory (DFT), a material where the layers are 
relatively strongly bound, at the T-point for up to five 
layers. 38,39 Using the scalar product of the displacement 


vectors as a metric for the similarity of the normal modes, 
we find a perfect match for the low-frequency modes, 
as those modes originate from acoustic modes. For the 
higher-energy optical modes, we find a small mixing ef¬ 
fect of modes with the same symmetry, which increases 
with the layer number. However, for five-layer WS 2 , we 
still find an average similarity of 93.4% of the linear-chain 
model normal modes compared to the DFT results. 

We expect all other layered materials to behave simi¬ 
larly; therefore, unless the displacement patterns are re¬ 
quired with high accuracy, the linear-chain model might 
be used to find the vibration patterns of the high-energy 
normal modes. 


IV. RESONANCE EFFECTS OF RAMAN 
MODES IN MOS 2 

We prepared single and few-layer M 0 S 2 flakes by me¬ 
chanical exfoliation on Si substrates with an Si0 2 layer 
of 50nm and 90nm thickness from natural M 0 S 2 (SPI 
supplies) and freestanding bilayer M 0 S 2 as described in 
Ref. 11. We determined the layer number by measur¬ 
ing the low-frequency Raman modes. 26,28,30 Raman mea¬ 
surements were performed with a Horiba LabRAM HR 
with 458 nm, 532 nm and 830 nm diode lasers, a T64000 
triple monochromator system with an Ar-ion laser and a 
HeCd-laser, and a Dilor XY800 spectrometer with a Kr- 
ion laser and a dye laser. The Raman measurements were 
performed at room temperature in a confocal setup using 
a lOOx objective except for the UV measurements, where 
a 40x objective was used. For all measurements the laser 
power was below 0.5 mW to avoid heating effects. 

As shown in Fig. 1, the new Raman modes in few-layer 
M 0 S 2 are observed at 286 cm -1 and at 471 cm -1 , in excel¬ 
lent agreement with our group-theory analysis (Table II). 
By polarization-dependent measurements in (x,x) and 
(x,y) scattering geometry, we confirmed our assignment 
of the 286 cm -1 peak to E g /E' symmetry and of the 
471 cm -1 peak to A\ g /A 1 symmetry. The spectra at 
2.7eV and 3.8 eV excitation energy shown in Fig. 1 (a) 
and (b), respectively, indicate different resonance behav¬ 
ior of the different Raman modes. Therefore, we compare 
in Figs. 8 and 9 the intensity behavior of the always ob¬ 
served Raman modes with that of the newly observed 
modes. 

Figure 8 shows the intensity of the Ai g /A 1 mode 
(407 cm -1 ) in single and bilayer MoS 2 , normalized to the 
number of layers and to the Raman intensity of the Si 
521 cnC 1 peak of the underlying substrate, as a func¬ 
tion of excitation energy. As expected, the M 0 S 2 Raman 
signal is resonantly enhanced at the A and B excitonic 
transitions. 40 Above 2.4 eV excitation energy, the Raman 
intensity increases even further, which we attribute to 
resonance with the optical absorption around the C ex- 
citon transition. This optical absorption peak has been 
recently determined to be at 2.84 eV in single-layer M 0 S 2 
and at 2.73 eV in bilayer M 0 S 2 41 ; the intensity increase 







Excitation energy (eV) 


FIG. 8: Intensity ratio of the normally observed 
Raman modes and the 521 cm” 1 mode of the silicon 
substrate, normalized to the number of M 0 S 2 layers. 
Hollow (filled) symbols indicate the Alg/A^ « 

407 cm -1 (E 2 g /E' « 383 cm” 1 ) Raman mode. A, B , 
and C denote the excitonic transitions. Different 
thickness of the SiC >2 (50 nm or 90 nm) are indicated. 
For the freestanding bilayer sample the Raman intensity 
of the surrounding substrate was used to calculate the 
intensity ratio. The lines connecting the data points are 
a guide to the eye. Data points in the range of the A 
and B excitons are taken from Ref. 40 


towards 2.7-2.8eV, as shown in Fig. 8, is in good agree¬ 
ment with these values. The second typically observed 
Raman mode, the E g /E' mode at 383 cm” 1 , shows a very 
similar intensity dependence on the excitation energy, as 
shown for 1L in Fig. 8. In order to check for possible 
effects of the substrate in the excitation-energy depen¬ 
dence of the Raman intensity, e.g., due to interference 
in the SiC >2 layer, we measured 1L and 2L M 0 S 2 on Si 
substrates with SiC >2 thicknesses of 50 nm and 90 nm, as 
well as freestanding 2L MoS 2 - Although the intensity ra¬ 
tios to the Si Raman intensity vary slightly, the overall 
intensity increase towards the C exciton is similar in all 
samples. 

The newly observed modes, on the other hand, start 
to appear in the spectra only above 2.4 eV excitation en¬ 
ergy. Their relative enhancement towards the C absorp¬ 
tion peak, however, is even stronger than that of the other 
Raman modes. This is shown in Fig. 9, where the inten¬ 
sity ratio between the “new” E g /E' mode at 286 cm” 1 
and the peak at 383 cm” 1 of the same symmetry is de¬ 
picted for 2L, 3L, and 4L M 0 S 2 . In 2L M 0 S 2 , this rela¬ 
tive intensity shows a maximum at ~ 2.7eV; in 3L and 
4L M 0 S 2 this maximum gradually shifts to lower excita¬ 
tion energy. This fits again nicely with the layer-number 
dependence of the C absorption peak 41 and supports our 
assignment to this resonance. At 3.8 eV excitation energy, 
the “new” Raman modes show even higher intensity than 



FIG. 9: Intensity ratio of the newly observed E’/E\ g 
Raman modes of M 0 S 2 and the E'/E\ g at « 408 cm” 1 
The lines connecting the data points are a guide to the 
eye. 


the other, typically observed modes, see also Fig. 1. The 
other newly observed mode (A- tg /A\ at 471cm” 1 ) shows 
a similar trend. However, it overlaps with the second- 
order Raman spectrum of the LA and LA' phonons from 
the M point 42 , which makes an analysis of the intensity 
less exact. In monolayer M 0 S 2 the new modes are not 
observed, independent of excitation energy, and in agree¬ 
ment with the symmetry considerations discussed above. 

In order to understand the distinct resonance behavior 
of the newly observed Raman modes, we consider the 
spatial distribution of the excitonic wave functions in few- 
layer M 0 S 2 . Theoretical calculations and experimental 
results predict the wave function of the A exciton in few- 
layer M 0 S 2 to be strongly confined to a single layer with 
only small overlap to the neighboring layers. 43,44 Recent 
calculations of few-layer MoSe 2 , which is quite similar in 
electronic structure to M 0 S 2 , also showed that the wave 
functions of the conduction and valence band states at 
the K point are strongly confined to a single layer. 45 

As we will discuss below, our results suggest that the 
C exciton, on the other hand, is much less confined to a 
single layer in few-layer M 0 S 2 , but instead extends over 
the entire few-layer thickness. This is in agreement with 
the significantly larger blueshift of the C absorption peak 
compared to the A and B absorption peaks when the 
number of layers is reduced 41 . 

The spatial confinement of the wave functions has 
strong implications on the Raman scattering process. If 
the Raman excitation is in resonance with the A or B 
exciton in few-layer M 0 S 2 the excitation is still confined 
mainly to a single layer [see Fig. 10(a)]. As a result, the 
A-layer system reacts rather like an independent super¬ 
position of N single layers from a selection-rule point of 
view . In the single layer, the vibrations at 286 cm” 1 
and 471cm” 1 are symmetry-forbidden in backscattering 
geometry, and therefore their corresponding symmetry- 
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(a) A lg - Raman active 



(b) A 2u - Raman inactive 





A 2 "- Raman inactive 



A,'- Raman active 



. M 


FIG. 10: In resonance with the A exciton, the Raman 
active A\ g mode of 2L M 0 S 2 can be “seen” like the 
Raman inactive A 2 mode of 1L M 0 S 2 (a). The Raman 
inactive A 2u mode of 2L M 0 S 2 can be “seen” like a A lg 
symmetry mode of 1L M 0 S 2 (b) and starts to become 
observable. 



FIG. 11: Schematic views of the optical absorption 
process (black arrow) for (a) 1.5 eV, (b) 1.9 eV, (c) 

2.8 eV and (d) 3.8 eV excitation energy. The red lines 
show the difference in energy of the conduction band 
(E c ) and valence band (E„) states, taking the spin-orbit 
splitting of the conduction band into account. The 
black lines indicate the exciton ground states. The gray 
dashed line indicates a virtual transition. 


allowed Raman signal in the few-layer system is still van¬ 
ishingly small in resonance with the A or B exciton. Only 
when the optical excitation creates an exciton that is ex¬ 
tended over the entire number of layers, as in resonance 
with the C transition, the symmetry of the few-layer 
system becomes dominant and determines the selection 
rules, allowing the newly observed Raman modes and 
bringing them into resonance. 

The different regimes of optical excitations are 
schematically illustrated in Fig. 11. We show the differ¬ 


ence between valence-band and conduction-band energies 
schematically from a density-functional theory (DFT) 
calculation of the single-layer M 0 S 2 band structure. Be¬ 
cause the DFT results underestimate the band gap, we 
applied a constant scaling factor to the energies such that 
the A exciton is at 1.9 eV when assuming a binding energy 
of 445meV. 46 Figure 11 (a) illustrates a non-resonant ex¬ 
citation at 1.5 eV, which leads to very weak Raman signal. 
We found it to be even weaker than the Si second-order 
modes around 300cm -1 (not shown here). 

At 1.9eV [Fig. 11(b)], the A exciton resonance is 
reached; because of the intralayer confinement, the Ra¬ 
man spectra of few layers and single layers are rather 
similar. Note, however, that the van-der-Waals and elec¬ 
tronic interactions between the layers do exist and lead 
to the observed shifts in the Raman frequencies. 15,4 ' At 
« 2.8 eV [Fig. 11(c)], the C exciton is reached, which ap¬ 
pears to be the first resonance which is extended over all 
layers. In this case, the full symmetry of the few-layer 
system is “seen”, and all symmetry-allowed vibrational 
modes become resonant (Table II), in particular the ones 
that are not Raman active in the single layer or bulk, i.e. 
the new modes at 286 cm -1 and 471cm -1 . 

Above this energy, we observe the same situation, i.e., 
the optical excitation is not any more confined to individ¬ 
ual layers. Figure 11(d) indicates a possible resonance for 
3.8 eV excitation energy, which might be in the contin¬ 
uum of optical excitation at the nearly flat bands between 
the T and K point of the Brillouin zone. 

Our results suggest that a high intensity of those Ra¬ 
man modes that are specifically only allowed in few-layer 
systems (compared to their single-layer or bulk counter¬ 
part) in general indicates an optical excitation that is not 
confined to a single layer in the few-layer system. Vice 
versa, the absence of these modes indicates ( i ) a true 
single-layer sample or (ii) a few-layer sample where only 
intralayer-confined states have been excited. 

Examples supporting our interpretation can be found 
in literature: In few-layer WSe 2 , where the A exci¬ 
ton energy is 1.6-1.7 eV, the new Ax g /A! 1 mode has 
been reported for excitation energies larger 2.33 eV. 9,18,19 
At lower excitation energies, the new mode was not 
detectable. 18 In few-layer MoTe 2 , the new Ai g /A’ l mode 
has been reported for excitation energies of 1.96 eV and 
2.33 eV, again in both cases the excitation energy are 
much larger than the optical band gap of <1.1 eV. 17,48 

Finally, resonant excitation into intralayer-confined 
excitons in few-layer systems can by itself “activate” 
Raman-inactive modes [see Fig. 10(b)]: As discussed 
above, then the selection rules follow more the rules of N 
independent layers. Thus, N Raman peaks correspond¬ 
ing to a given single-layer Raman mode can be observed 
- given that the frequency splitting between them can 
be spectrally resolved -, although only every other mode 
of the set of N modes is Raman active in the TV-layer 
system (Table II). In this case, the truly Raman-allowed 
modes still dominate the spectra; the A-layer inactive 
but lL-active modes are significantly weaker in intensity. 
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Such Raman spectra have been recently observed in WS 2 , 
which exhibits particularly large frequency splitting for 
the optical A-\ g jA\ modes. 39 

V. CONCLUSION 

In conclusion, we have presented a generalized 
symmetry-based analysis of vibrational modes in TV-layer 
materials. We have shown that the vibrational spectra in 
few-layer systems can have distinct signatures for the in¬ 
teraction between the layers compared to the correspond¬ 
ing single-layer or bulk material. Moreover, we have 
shown how to derive the vibrational modes, their sym¬ 
metry, displacement vectors, and approximate frequency 
range from the symmetry of the single-layer system. 

In general, two symmetry-related effects determine the 
vibrational spectra of few-layer systems: (i) the symme¬ 
try group of the few-layer system can be different from 
the single-layer or bulk symmetry group. This can lead 
to activation of vibrational modes, such as the silent bulk 
B- 2 g mode in the D§h symmetry group of M 0 S 2 , which be¬ 
comes a Raman active A lff mode in N layers with even N 
(symmetry group D^d)- (**) each vibration of the single 
layer will always lead to several combination modes of 
the few-layer system that show opposite behavior under 
inversion (horizontal reflection). As a result the combi¬ 
nation modes resulting from Raman inactive vibrations 


of the single layer can be combined into Raman active 
vibrations of the few-layer system, such as the symmet¬ 
ric breathing-like vibrations arising from Raman-inactive 
acoustic modes, or the optical A” mode in single-layer 
M 0 S 2 which becomes the Raman-active A\ mode in N 
layers with odd N. Our experimental results of such 
additional Raman modes in few-layer M 0 S 2 confirm our 
theoretical predictions. 

Based on the strong resonance of the few-layer specific 
Raman modes with the C optical transition, we predict 
that the corresponding exciton wave function is extended 
over all layers of the few-layer M 0 S 2 , in contrast to the 
A and B excitons, which appear strongly confined to a 
single layer. Consequently, the intensity of “new” few- 
layer Raman modes in layered materials in general can 
indicate the spatial distribution of the optical excitation. 
The symmetry-based analysis and the resulting predic¬ 
tions can in general be applied to any layered crystal 
structure with inversion symmetry and/or horizontal re¬ 
flection symmetry. 
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